Anisotropic electron gr-factor in quantum dots with spin-orbit interaction 
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g-factor tuning of electrons in quantum dots is studied as function of in-plane and perpendicular 
magnetic fields for different confinements. Rashba and Dresselhaus effects are considered, and 
comparison is made between wide- and narrow-gap materials. The interplay between magnetic 
fields and intrinsic spin-orbit coupling is analyzed, with two distinct phases found in the spectrum 
for GaAs in perpendicular field. The anisotropy of the g-factor is reported, and good agreement 
with available experimental findings is obtained. 
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Spin properties in semiconductor quantum dots (QDs) 
have become a field of intense research because of the 
possible use of the electron spin degree of freedom as a 
quantum biti It is then essential to have a clear under- 
standing of the processes that may induce spin relaxation 
on electrons in the QD; long spin relaxation times, and as 
pure as possible spin states, are required so that spin can 
indeed transport information without losses. Insights on 
the purity of the spin degree of freedom of electrons in 
QDs can also be extracted from measurements of their 
effective g-f actor. 

Manipulation of g-factor in semiconductors and de- 
scription of its tensorial nature have been considered 
theoretically as well as experimentally. Among vari- 
ous techniques, appropriate system design has achieved 
gate-voltage control of 5-factor of electrons in a quan- 
tum well^ and spin manipulation using gigahertz elec- 
tric fields."^ g-factor measurements have been reported 
in QDs by means of capacitance^ and energji^ spectro- 
scopies, for example. The tensorial nature of electron 
(/-factor in spherical QDs,^ as well as surface^, and spa- 
tial confinemenliS effects, where spin-orbit (SO) coupling 
plays an important role, have been addressed theoret- 
ically. We have recently reported on the influence of 
SO coupling on the electronic spectrum of 2D parabolic 
QDs.^" Spin-orbit effects on factor have been addressed 
for 2D electrons in a quantum welUi and for electrons in a 
parabolic QD,^^ with emphasis on the difference between 
narrow- and wide-gap compounds. 

One of the main causes of spin relaxation and g-factor 
variation is the SO interaction. When QDs are built in 
semiconductors of zincblende structure in the plane of 
a 2D system, there are two possible forms of SO cou- 
pling, namely the Rashba-'^'^ and Dresselhaus^'* interac- 
tions; the former is due to the surface inversion asym- 
metry (SI A) induced by the 2D confinement, while the 
latter is caused by the bulk inversion asymmetry (BIA) 
intrinsic in zincblende structures. The SO coupling mixes 
spins with different orientations in the Zeeman sublevels, 
which yields an intrinsic spin relaxation source and pro- 
duces variations in the QD g-factor from the pure Zeeman 
splitting expected in a magnetic field. 

In this work we study the anisotropy of effective g- 



factors in 2D parabolic QDs, and analyze the intricate 
competitioiii^ between external magnetic fields and in- 
trinsic SO couplings. Wide- and narrow-gap materials 
with different confinement potentials are considered un- 
der in-plane and perpendicular magnetic fields. We show 
that the f;-factor can be tuned to have positive, zero, 
and negative values at given perpendicular fields, and 
find that two distinct phases can be present in the QD 
spectrum, where at low (high) fields the SO coupling en- 
hances (suppresses) the Zeeman sublevel splitting. We 
find that the widely used perturbative approach-- based 
on a unitary transformationi^ of the SO Hamiltonian has 
strict limitations in QDs if the well defining the 2D con- 
finement has small width. We also show that even for 
in-plane magnetic fields - typically assumed to yield the 
bulk ^'-factor without orbital contributions - the SO cou- 
pling reduces the value of g, especially for large QDs. 

The QD is defined by an in-plane parabolic confine- 
ment, V{p) — mtjQ/9^/2, where m (wo) is the electronic 
effective mass (confinement frequency), while the per- 
pendicular confinement V{z) is assumed strong enough 
to reduce the electronic states to the first conduction 
subband of the quantum well; its function is ip{z) — 
■y/2/zo sin (ttz /zq) , zq being the QD vertical width for 
hard wall potential. In the absence of SO interac- 
tions, the QD Hamiltonian is Hq — h^\s? /2m -I- V{p) -\- 
gotJ-B^± • o'/2, where k — — iV + eA/(hc) is the canon- 
ical momentum, go is the material bulk g-factor, a 
stands for the Pauli matrices, /is is the Bohr magne- 
ton, and A — B±p{— sm(j),cos(j),0)/2 describes a per- 
pendicular magnetic field Bx = i3_L(0,0, 1), which lifts 
both orbital and spin degeneracies of levels. Hq yields 
the Fock-Darwin (FD) spectrumjiS Enia, (2n + \l\ + 
l)hn -t- lhujc/2 + gf)iiBB±az/2, with effective (cyclotron) 
frequency Q — y^uiQ -\- {ojc = eB±/mc), where 

n = 0,1,2,... and / = 0,±1,±2,... are respectively 
the radial and azimuthal quantum numbers. The FD 
eigenfunctions are 'i/ni{x,(p) — i?„|;[ (x)e*''^/\/27r, where 

in terms 



of associated Laguerre polynomials L„ , A = y/h/{mfl) 
is the effective QD lateral length, and x — p/\. 

The SO Hamiltonian added to Hq is Hso — Hsia + 
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Hbia- The SIA term, due to the full confining 
potential V{t) — V{p) + V{z) and with coupHng 
parameter a, is given by Hsia — ao" ■ V\^(r) x 
k. It can be separated as Hsia = Hr + Hlat, 
where Hr = —a/ X{dV/dz)[a+L^A^ + cr^L+A^] is 
the Rashba terrain caused by the interfacial electric 
field generated by the perpendicular confinement, and 
Hlat = a{fuvo/l^)a,[L^ + {X/lefx^ /2\ is due to 
the lateral confinement (diagonal in the FD basis); 
we define L± — e^"^, a± — [a^ ± z(Ty)/2, = 
-id/d(j>^ and A± = ^d/dx + L^/x + {\/IbYx/2, and 
the confining (magnetic) length is Iq = ^ hjiraujo) 
{Ib — \/fi/ (tooJc))- The z-confinement yields the BIA 
term^ Hbia = iWxkx {k^ - kl) + ayky {kl - kl) + 
Czkz (fc^ — ky^] = Hj^ + Hj^\ the linear contribution is 
H^ = -i^ (kl) /\[cF+L+A+ -a^L^AJ\, while the cu- 
bic is H% = i-i/\^[a^L\Hi + cr+L?_7?2 + a.L.H^ + 
a+L+Hi], where H^ ^ A^ + [X/IbTB, + (A/Zb)^^ + 
(X/Ib)^ Di, with i = 1,2,3,4. The sixteen functions Ai, 
Bi, Ci, Di are known^ 7 is the coupling parameter, and 
(^z) = (""/^o)^ for hard wall confinement. The total QD 
Hamiltonian, H ^ Hq + Hlat + Hr + H^ + is 
diagonalized in a basis set including 110 FD states. 

When considering an in-plane field (i3||) we take the 
full Hamiltonian above at zero i?_L-field, H{B± = 0), 
plus the Zeeman contribution go^J-B'B\\ • f/Z, which lifts 
only the spin degeneracy of levels. 

Before discussing the QD effective factor, we com- 
ment on the infiuence of SO coupling on the spectrum 
and its competition with external magnetic fields. This is 
done in Figs. ^ and [3 respectively for InSb^^ and GaAs^ 
QDs, the former having larger SO coupling and Zeeman 
splitting than the latter. In these two figures, left (right) 
panels refer to an in-plane (perpendicular) field, and dot- 
ted lines in panels A and B refer to the pure FD spectrum 
(no SO-coupling). In panels A and B of Fig.^ the main 
effects of the SO interaction on the InSb QD levels are as 
follows: crossings between FD levels are converted into 
anticrossings (AGs) according to the selection rulesiSi of 
Hso (and shifted to higher field values in the perpendicu- 
lar case) ; the low-field spectrum is displaced to lower en- 
ergies, while at high-field is not altered much; zero-field 
splittings appear in the spectrum, and the original FD 
energy shells are separated into components according to 
the total angular momentum projection, j — / + crz/2. 
Panels C and D (E and F) of Fig. □ show the field- 
evolution of spin (az) (orbital (l)) angular momentum 
expectation values for the seven lowest QD levels, where 
it is clearly visible how the SO coupling mixes states, re- 
sulting in (Tz and I no longer refiecting pure FD levels, 
especially around AC points. For the purposes of dis- 
cussion, we still label states by the numbers {l,az} even 
under SO interaction. 

By comparing left (right) panels in Fig.^ one can ob- 
tain the QD level sequence and realize which SO mecha- 
nism is responsible for the level AGs under (B±). Fo- 
cusing on the lowest energy levels, notice at ~ 0.1 T 



FIG. 1: InSb QD spectrum [isll . Left (right) panels refer 
to in-plane (perpendicular) magnetic field. Panels A and B 
show spectra (dotted lines refer to QD without SO interac- 
tion). Panels C and D {E and F) show field-evolution of spin 
(orbital) angular momentum expectation values for the seven 
lowest levels; ground state is spin-up. Labels in C and D 
show which SO term is responsible for the lowest AGs. Col- 
ors indicate different levels: black, red, green, ... stand for 
ground, first excited, second excited, ... states, in all panels. 
A sudden change in line color indicates a level crossing. 

FIG. 2: Level spectrum for GaAs QD |20ll . Bx-scale of panel 
B is different from panels D and F. No AG occurs under By . 
Ground state at low _Bx-fields is spin-down, and it is flipped 
back to spin-up around 13.7 T, as shown in panel D. 



that the lowest two states in the first shell have \i\ — 1/2 
{{l^Gz} = {0,1}, {0,-1}); in the second shell, the two 
lower states have \j\ = 1/2 ({—1,1}, {1,-1}), while the 
two upper ones have \j\ = 3/2 ({1,1}, {—1,-1}). At 
B±^ ~ 0.1 T, the only difference in such sequence is in 
the ordering of the two upper states in the second shell, 
which become \j\ = 3/2 ({-1,-1}, {1,1}). This order- 
ing reversal is due to the influence of B±^ on the orbital 
features of the state, which changes the effect of SO cou- 
pling on these levels; similar features appear in higher 
shell states. Even though the ground state does not ex- 
hibit AG for any field-direction, it is also not spin-pure at 
low magnetic fields; observe that (ctz) ~ 0.75 in both B± 
and B\\ , and it is positive since go is negative. Regarding 
level AGs, the lowest under (~ 4.8 T) involves the 
second and fourth QD states in panel A, which at that 
field are the levels {0,-1} and {1,1}; such levels are con- 
nected by the operators cr±L±, and although they appear 
in both H'j^i and Hj^, terms, this AG is mostly due to the 
linear Dresselhaus term as indicated in panel C . On the 
other hand, the lowest AG under B± (~ 3.3 T) involves 
the second and third QD states in panel S, which at 
that field are the levels {0,-1} and {—1, 1}; such levels 
are connected by the operators cr±Lip, so that this AG 
is due to the Rashba term Hr, as indicated in panel D. 
Higher energy AGs are evident at similar values of -By 
(or B±) and are due to the same SO mechanism. The 
different SO terms in H, as well as the level dispersion 
in the two field-configurations producing the AGs, result 
also in different spin mixings. This clearly affects the QD 
effective g-factor, as we will see below. AGs due to the 
cubic Dresselhaus term are also observed at higher fields: 
at B± ~ 6 T the lowest AG involves the fourth and fifth 
levels, {0,-1} and {—3,1}, connected by the operators 
cr^Lj.. Interestingly, at ~ 14 T, where AGs due to 
H]^ involve higher energy levels, (/) for all QD levels col- 
lapse to zero, indicating strong orbital mixing with full 
spin polarization. 

Figure 121 has the same analysis for a GaAs QD. The 
left panels show that in an in-plane field the SO cou- 
pling is not strong enough to induce AGs in the spectrum 
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(panel A), even though zero- field splittings are clearly ob- 
served for different |j| values (e.g., look at the second 
shell). Notice that levels present the same ordering as 
in the low-i3|| InSb QD, and that spin mixing is present 
in the spcctrmn, as states do not have (az) — ±1 (panel 
C); e.g., the ground state has (ctz) ~ 0.75 at small fields. 
This shows the intricate interplay between external fields 
and SO coupling in the definition of electronic properties 
of QDs:^^ In a perpendicular field (right panels of Fig. 
O, however, a new feature is observed: the SO interac- 
tion flips the spins of the two lowest levels at low fields, 
so that the QD ground state becomes a spin-down level. 
For this reason, the lowest AC at B± ~ 2 T (notice differ- 
ent B^-scale in panel B compared to D and F) involving 
the second and third levels, {0, 1} and {—1,-1}, is due 
to the linear Dresselhaus term. Observe in panel D that 
those two lowest levels have the same value of (cr^) ~ —1 
between 4 and 8 T, and that the ground state flips back to 
its expected spin-up character of a go < QD around 13.7 
T, when a level crossing is verified in the spectrum (not 
shown) . This result prompts one to consider two distinct 
ways of defining the QD effective g-factor, and points out 
strict limits for the validity of the widely used perturba- 
tive approach^ (or unitary transformationii) for deal- 
ing with SO effects in QDs, especially if the well thick- 
ness zq defining the z-confinement^-'' is small. One has 
to keep in mind that these results depend on the QD en- 
ergy (length) scales: higher Eq — huo (smaller ^o) QDs 
are expected to have weaker SO effects, since their level 
spacing becomes larger. 

What is the 'correct' way of defining the effective elec- 
tron ^-factor in a QD? There are two possible defini- 
tions involving the two lowest Zeeman sublevels for each 
field-direction, namely, g^^^/go = ^E/{goiJ,BB±^\\) or 
9±.\\/9n = (Acrz)/2, where AE ((Acr^)) is the energy 
splitting (spin expectation value difference) of those lev- 
els under Although the first definition is used op- 
erationally in experiments where AE is measured, the 
latter is intuitively reasonable since the ^-factor is a 
quantity intrinsically related to the spin value of those 
levels. For no SO interaction, both definitions yield 
9± — 9\\ — ga, so that no anisotropy is present in the g- 
factor (other than the material anisotropies)i2i4 Figures^ 
(InSb) and0](GaAs) present the factor for different QD 
lateral sizes Zq, with left (right) panel for in-plane (per- 
pendicular) fields, and dotted lines show results without 
SO coupling. Panels A and B [C and D) use the defini- 
tion of g in terms of (Aa^) (AE), whose values can be 
inferred from panels C and D of Figs. ^ and |3 (panels E 
and F of Figs. El and Ej). Notice that for InSb in Fig. 
both definitions yield basically the same low-field results 
for the case (panels A and C). The drop in is faster 
in the g'^' curves because the two lowest states have the 
same spin at higher fields, while AE reaches a constant 
value (panel E); dotted lines show the field value where 
the original FD lowest crossing is converted into an AC 
by the SO coupling at a given Eq- 

In the low- field B± case (panels B and D), however. 



FIG. 3: Electron jr-factor of InSb QDs with different confine- 
ment energies Eq (or lateral sizes lo). Left (right) panels refer 
to in-plane (perpendicular) field. Panels A and C {B and D) 
show {g±), while panels E and F show sublevel splitting 
AE for the two lowest QD states. Panels A and B (C and D) 
show g as obtained from the (Actz) (Ai?) calculation; values 
of (Acr^) for the 15 meV QD are obtained from Fig. Q Dot- 
ted lines in all panels refer to corresponding QDs without SO 
coupling. AE in panel F is multiplied by 2. Color legend is 
shown in panel B; dotted lines follow same scheme. 



FIG. 4: Same as Fig. El but for GaAs QDs. g"^' has inverted 
sign at low fields indicating that SO coupling is stronger than 
the Zeeman splitting ((Acr^) of the 1.1 meV QD from Fig. 
EJ. g± has large values at low fields. Dashed lines connecting 
panels D and F show fields Bg for which g±/go = 1, and the 
SO effect is effectively cancelled. Color legend is in panel A. 
g in panel D is divided by 10. 

the two g-factor definitions yield different values mainly 
at weaker confinements, although the drop in g± is also 
faster in the g"^^ curves {AE values are found in panel 
F). Observe that for the smallest confinement energy 
(largest SO effect) of 3.0 meV, a sign change is seen in 
g'^' around 1 T, which relates to an unusual crossing in- 
volving the ground and first excited states (compare with 
panel B of Fig. ^ where this crossing is absent). In both 
field-directions, smaller Eq yields smaller g-factor, which 
shows that SO coupling provides a channel to manipulate 
g in QDs under magnetic fields. At the same time, a mea- 
surement of effective g under i?_L.|| might give information 
about the values of SO constants since, as discussed in 
Fig. n the lowest AC in the spectrum of QDs in differ- 
ent field-directions is due to a distinct SO mechanism. 
Notice the clear anisotropy in (panels C and D): the 
same QD confinement shows gf^ < g^ , since the mixing 
with higher orbitals is stronger for B± . If the g"^' defini- 
tion is considered (panels A and B), such anisotropy is 
not as remarkable at low fields, although a sign change 
is obtained at E'o = 3 meV, only in the g"^ curve. No- 
tice that unlike the case with no SO interaction, even 
-B|| reduces the ^-factor and this reduction can in fact be 
strong (> 50 % at Iq ^ 300 A). 

FigureElshows that for GaAs both definitions in panels 
A and C give essentially the same values of , where AE 
is shown in panel E. Differences are noticed only at high 
magnetic fields for the weakest confinements (0.7 and 1.1 
meV). Results are totally different for the perpendicular 
field case. Panel B can be understood by looking at 
panel D of Fig. |2 (for 1.1 meV), where it can be seen 
that g"^' has inverted sign at low fields, becomes zero for 
B± between 4 and 8 T, then acquires inverted sign again, 
and suddenly flips back to its 'normal' behavior at 13.7 
T; under even higher fields, g"^' goes to zero since the two 
lowest level spins are aligned. Panel B also shows that 
larger Eq values (thus smaller SO coupling) cancels the 
field range where g"^' is zero. Notice that the weaker the 
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confinement the smaller the field where the sign change 
occurs. 

One finds totally different results for gf^ (panel D, with 
AE shown in panel F), which may even assume values 
12 times larger than 170 for the smallest Eq at low fields; 
still at low fields, larger Eq tend to reduce towards 
go. At high fields (inset G), gf goes to zero when the 
level crossing involving the ground state occurs, and af- 
ter such field it increases again. For every Eq there is a 
magnetic field Bg - indicated by the dashed lines con- 
necting panels D and F - where gf^ goes from higher to 
smaller values than g^. In Panel F, differently from what 
occurs in panel E (and in panels E and F of Fig. ISJ, the 
Bg-field defines two distinct phases in the spectrum for 
a given Eq: below (above) Bg the SO coupling increases 
(decreases) the Zeeman sublevel splitting as compared to 
the case without SO interaction; one can then say that 
at B± — Bg, the SO coupling is cancelled by the mag- 
netic field in the formation of the sublevel splitting. Such 
result emphasizes the intricate competition between ex- 
ternal magnetic field and intrinsic SO coupling in QDs. 
[Notice that if a broad Gaussian well or a larger zq is con- 
sidered, such phases are not observed i^-''] In GaAs QDs, 
the anisotropic nature of the (7-factor is much more pro- 
nounced, despite the small values of the SO constants. 

We conclude with an experimental comparison. Ref. j^l 



reported Zeeman sublevel splittings in GaAs QD having 
Eq = 1.1 meV under in-plane field. They found AE ~ 
200 fieV at 10 T, while the corresponding curve in panel 
E of our Fig. ^yields AE ~ 180 ^eV. In a linear fit, they 
found \g\ — 0.29 ± 0.01, and from panel A (panel C) of 
Fig. Hone has |.gf[-'| = 0.30 (|.gf| = 0.31) at S = 10 T, 
while |(7[| I — 0.30 is found a,t B — from both definitions. 

We have shown how SO coupling is able to tune the 
electron ^-factor in QDs and even change its sign. We 
have analyzed the interplay between SO and Zeeman 
splittings on QD spectra and shown which SO term 
causes the lowest AGs in in-plane and perpendicular 
fields, as well as in different materials. We have seen 
that for GaAs QDs under B± the ground state has its 
spin character inverted at low fields if a narrow well con- 
fines the system in the z-direction. We have identified 
phases of B± in the spectrum where the SO interaction 
increases or decreases the Zeeman splitting of the lowest 
QD levels, and explicitly shown the anisotropic nature 
of (7- factor in QDs. All these features would not have 
been accessed if a perturbative approach had been used, 
especially for QDs with large lateral size. 

We thank support from NSF-IMG grant 0336431, 
GMSS at OU, and the 21"* Gentury Indiana Fund. 
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